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We investigate the thermalization via heat radiation of a silica fiber with a diameter smaller than 
the thermal wavelength. The temperature change of the subwavelength-diameter fiber is determined 
through a measurement of its optical path length in conjunction with an ab initio thermodynamic 
model of the fiber structure. The results differ significantly from the predictions of Planck's law 
based on the spectral emissivity of silica. Excellent agreement is obtained with a theoretical model 
that considers heat radiation as a volumetric effect and that that takes the emitter shape and size 
relative to the emission wavelength into account. These results are of fundamental interest, may 
lead to technical applications, and can contribute to improved models of the earth's climate system. 



Thermalization via heat radiation is an omnipresent process which, e.g., determines the temperature of stars and 
planets or the functioning of incandescent lamps. For a perfectly black body, the spectral emissive power of thermal 
radiation was first explained by Planck [1 who used quantized energies for the radiation field, thereby breaking the 
grounds for quantum theory. However, a real object does not radiate like a perfectly black body but emits a modified 
spectrum which can be related to Planck's law by introducing a correction factor, the so-called spectral emissivity 
[2]. Usually, this spectral emissivity is treated as a surface property that depends on the specific material and the 
surface roughness. Yet, if the material is not opaque, thermal radiation from within the bulk material can escape into 
free space provided that its point of emission is closer to the surface than its absorption length in the material. As 
soon as the latter gets larger than the size of the radiating body, one therefore has to consider thermal radiation as a 
volumetric effect [3J H] • 

Recent experiments have shown that size and geometrical structure of an object may play an important role for 
the spectrum, coherence, and angular distribution of its thermal radiation [5-9_. Moreover, it has been shown that 
the radiative heat transfer rate can be strongly enhanced by near-field effects [10]. Here, we optically measure the 
thermalization of a silica nanofiber with a diameter of 500 nm, smaller than all absorption lengths of the thermal 
radiation spectrum. The fiber resides in an ultra high vacuum chamber and heat transport is dominated by far-field 
radiation. Our data shows, to the best of our knowledge for the first time, that the total far-field radiated power of a 
subwavelength structure deviates from Planck's law and agrees quantitatively with the predictions given by so-called 
fiuctuational electrodynamics (FED) 

The silica nanofibers used for our experiment are realized as the waist of a tapered optical fiber (TOF) which is 
produced from a standard optical fiber in a heat-and-pull process |12l H3| . A temperature change of the nanofiber 
waist induces a change of its optical path length which is primarily due to the thermo-optic effect of silica. In order 
to measure this optical path length change, the TOF is enclosed in a fiber-based Fabry-Perot-type optical resonator 
which uses two fiber Bragg-gratings (FBGs) as cavity mirrors, see inset of Fig. [I] |14| . In this configuration, the 
optical path length change is translated into a frequency shift of the Fabry-Perot resonances which we read out via 
the transmission of a probe laser field, see Fig. [T] The nanofiber waist is heated by sending a second laser field 
through the TOF. Its wavelength lies in the transmission band of the FBGs and its power is adjustable by means of 
an acousto-optic modulator (AOM). A small fraction of this optical power is absorbed along the nanofiber, thereby 



FIG. 1. Schematics of the experimental setup for measuring the thermally induced optical path length change of the optical 
nanofiber. The inset shows the TOF-resonator with the nanofiber waist. 
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FIG. 2. Thermally induced optical path length change of the silica nanofiber. Circles: Experimental data; blue bands: FED 
predictions; green bands: predictions using Planck's law. (a) Example of thermalization dynamics, (b-d) Maximum optical 
path length change as a function of the heating laser power transmitted through the TOF and of the background gas pressure, 
respectively. 



heating the latter. After exiting the TOF, the transmitted probe and heating laser fields are separated by a dichroic 
mirror (DM) and spectral filters such that their respective powers can be measured independently using photodiodes 
(PD1 & PD2). By abruptly switching the heating laser beam with a mechanical shutter, the heat source can be 
turned on and off, allowing us to measure the heating and cooling dynamics of the nanofiber. 

For this purpose, we record a time trace of the probe transmission which consists of a sequence of Fabry-Perot 
transmission peaks. From one peak to the next, the single-pass intra-cavity optical path length has changed by half 
the vacuum wavelength of the probe laser, AL opt = Ao/2 = 426 nm. Prior to each cycle of heating and successive 
cooling, we establish a reproducible thermal situation by an initializing pulse sequence of the heating laser: First the 
heating laser power is set to a high value (80-160 mW) and is sent through the TOF for 30 s. This heating time is 
significantly longer than the thermalization time constant of the unprocessed fiber parts of ^8 s, determined in an 
independent measurement by tracking the central wavelength shift of the FBG stop-band. Then, the heating laser is 
switched off for 2 s. This cooling time is ten times longer than the thermalization time constant of the nanofiber, see 
below. This ensures that the nanofiber is thermalized at ambient temperature while the unprocessed fiber exhibits a 
constant initial condition. 

A typical time trace of AL opt is shown in Fig. [2^i for a TOF-resonator with a waist diameter of (500±50) nm, a 
nominal homogeneous waist length of L W aist = 10 mm, and a finesse of 85 ± 1 (TOF-resonator #1: detailed radius 
profile in the supplementary informations). For this measurement, the heating laser is switched on between t = s 
and t = ls, the transmitted power is Phoat = 32.7(1) mW, and the background gas pressure is p — 10~ 6 mbar. For 
both the heating and for the cooling process, AL opt varies rapidly during the first 500 ms. Following this initial fast 
dynamics, the variation of AL opt takes place on a much longer timescale on the order of several seconds. Because 
of the discrete nature of the measurement method which relies on counting transmission peaks, one additional data 
point was added at the end of the heating (cooling) period. Its value is chosen A /4 higher (lower) than the last 
measured data point and its error bar is given by ±Ao/4, thereby accounting for the fact that AL opt changed by less 
than Ao/2. 

In order to gain a quantitative understanding of the observed dynamics of AL opt , we develop an ab initio thermo- 
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dynamical model of the TOF including heat transport via thermal radiation, heat conduction along the TOF, and 
heat diffusion through the surrounding background gas (see supplementary informations) . The parameters that enter 
into the model are the specific heat, the heat conductivity, and the density of silica as well as the molecular mass and 
the number of degrees of freedom of molecular nitrogen in order to account for the background gas. Furthermore, 
the model uses the radius profile and the total radiated power per unit length of the TOF. The latter can either be 
computed according to Planck's law including the spectral emissivity of silica or using FED |15[ [16] . It has been 
shown that the transmission loss of TOFs is primarily governed by surface pollution [17j . In our model, we therefore 
assume that the heating of the TOF is caused by surface absorption, meaning that the heating power per unit length 
at any position along the TOF is proportional to the radius-dependent surface intensity of the TOF-guided mode 
|18| . integrated over the circumference. Furthermore, we neglect temperature gradients in the background gas. In our 
system, this assumption is valid for pressures p < 10~ 3 mbar. For higher pressures, our model may overestimate the 
effect of heat diffusion. In the low pressure regime, the largest uncertainties that enter into the model are the radius 
profile, which is known with a relative error of ±10 % [191120) . and the wavelength-dependent complex refractive index 
of fused silica, n. The latter is the sole material property that is used in the FED calculation of the radiated power as 
well as in the calculation of the spectral emissivity. It was extracted from various literature sources for wavelengths 
ranging from ^300 nm to ^50 fim (see supplementary informations). 

We numerically solve the two resulting differential equations (Planck and FED) for the same heating and cooling 
cycle as in the experiment and determine the time-dependent temperature profile along the TOF as a function of the 
total absorbed optical power, P a b s , and the background gas pressure. Four different parameter sets are used for these 
calculations, combining the minimum or maximum TOF radius profile with the values of h that yield the minimum or 
maximum radiated power per unit length, respectively. From these temperature profiles, we then compute AL opt via 
the temperature-dependent effective refractive index which is integrated along the TOF-resonator (see supplementary 
informations). The theory bands in Fig. are delimited by the two extremal time traces of AL opt that reproduce 
the maximum observed optical path length change AL™^ X (solid lines) . We find excellent agreement between the time 
dependence of AL opt predicted by FED and the experimental data while the prediction for Planck's law significantly 
deviates from the latter. 

Figure |2j) and c show Ai™^ x as a function of the heating laser power transmitted through the TOF, Pheat, for 
TOF-resonators #1 and #2. The latter corresponds to a TOF with a waist diameter of (500±50) nm, L W aist = 5 mm, 
and a finesse of 18 ± 1. In both cases, AL™* increases roughly linearly with Pheat within the observed power range. 
In order to analyze the measured data, the ratio rj between Ph ea t and P a bs has to be determined. This is done by 
inserting P a b s = V Pheat into both theories and fitting their predictions independently to the data using the four 
parameter sets concerning the radius and the complex refractive index. We compute the mean of the minimum and 
maximum resulting 7/- value, fj — (?7max + ?7min)/2, for both theories and plot the extremal theory prediction (solid lines 
delimiting the colored bands) using the corresponding two parameter sets and P a b s = V Pheat- For both resonators, 
the FED prediction is in excellent agreement with the experimental data while the prediction for Planck's law exhibits 
a significant qualitative deviation. 

From the temperature profiles along the TOF which lead to the two extremal predictions of AL™^, we determine 
two corresponding extremal predictions for the temperature at the center of the waist, Tfji'-J. and T™^ st , respectively. 
The resulting non-linear temperature scale, T waist = (T™ t + T^J/2, where and AT waist = [T™ x - T™£ t )/2, 
is shown on the right axes of all panels of Fig. [2] with the ticks and the gray bars indicating T W aist and AT wa ; s t, 
respectively. It covers a temperature range from room temperature to almost 2000 K. While this upper temperature 
limit exceeds the glass transition temperature of silica of about 1450 K [ST] , the corresponding viscosity of 10 6 Pa • s is 
still very high and ensures the mechanical stability of the TOF. We note that the temperature scale in Fig. [2^ refers 
to the cooling process. The corresponding scale for the heating process differs from the latter and is not shown for 
clarity. 

Figure|2jl shows AL™^ X as a function of the background gas pressure for TOF-resonator #2 and Ph ca t = 6.7(2) mW. 
Both the measured values and the theory prediction yield constant results for pressures lower than 10~ 4 mbar, 
confirming our assumption that heat diffusion through the background gas is negligible for all other measurements 
presented here, which were performed at pressures below 10~ 6 mbar. 

We now analyze the thermalization dynamics for varying heating powers and thus for varying AL™^ 51 . For this 
purpose, we determine the 10 % - 50 % and 75 % - 90 % rise times (heating) and the 90 % - 50 % and 25 % - 10 % 
fall times (cooling) of the time traces of AL opt in order to quantify the initial and the final dynamical behavior. We 
find time constants on the order of 100 ms which are plotted in Fig. [3] as a function of AL^ X . The same quantities 
are derived from the theory predictions for Planck's law and FED, respectively, and plotted as bands in the same 
graphs. The agreement between the experimental data and the FED predictions is excellent, whereas Planck's law 
predicts up to one order of magnitude shorter time constants. All theory predictions in this figure are ab initio results 
without any adjustable parameters. 

Summarizing, we present conclusive experimental evidence for the failure of Planck's law in predicting the total 
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FIG. 3. Thermalization time constants of the silica nanofiber as a function of the maximum optical path length change for the 
two TOF resonators. Blue circles: measured initial time constants; red squares: final time constants; blue (red) bands: FED 
predictions for initial (final) time constants; green (yellow) bands: predictions for initial (final) time constants using Planck's 
law. 



thermally radiated power of a silica fiber with a sub-micron diameter. Moreover, we quantitatively confirm that 
its heat radiation is a volumetric effect and that fluctuational electrodynamics yields an accurate description of this 
phenomenon. These results have important implications for technical applications, such as heat management in 
nano-devices or the realization of energy-efficient incandescent lamps. Moreover, they might also lead to a better 
understanding of the impact of particulate matter, such as mineral dust aerosols from soil erosion and soot from 
combustion sources, on the climate system via absorption and emission of solar and thermal radiation [22] [33]. 
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SUPPLEMENTARY INFORMATIONS 



Optical path length change 

The optical path length change has its origin in three basic effects: The thermo-optic and the strain-optic effect of 
silica, i.e., the temperature- and strain-dependence of its refractive index, n, as well as the radial thermal expansion 
of the fiber. All effects lead to a change of the so-called V-parameter of the vacuum-clad fiber waist [18 , V — 
2-Ka\/n 2 — 1/Ao, either via a change of n or via a change of the fiber radius a. This results in a change of the effective 
refractive index of the fundamental HEn fiber mode in which the light is guided through the TOF, n e g = /3{V)/ko, 
where f3(V) is the V-dependent propagation constant of the HEn mode [181 and fco = 2tt/\o. 

The tensile stress of the fiber varies with temperature because of the axial thermal expansion of the fiber: Following 
the heat-and-pull process, the TOF is pre-strained before fixing it to its mechanical mount, leading to an elongation 
of the nanofiber waist by AL str ain/^o = 5 • 10~ 2 of its original physical length Lq. This elongation is much larger 
than the axial thermal expansion of the nanofiber which, assuming an unstrained nanofiber and a thermal expansion 
coefficient of silica of a tc = (AL/L) /AT = (4- 7) • 1CT 7 /K [24], would be of the order of AL theim /L « 1(T 3 for the 
maximum measured temperature of 2000 K. As a consequence, the axial thermal expansion of the nanofiber does not 
modify the length of the pre-strained nanofiber but rather lowers its tensile stress. This leads to a change of n via 
the strain-optic effect and to an increase of a due to the reduced transversal contraction of the nanofiber. 

The thermally-induced change of the effective refractive index for a given fiber radius is thus given by: 

dn eS (a,T) dn cS (a) ( dn(T) \ dn eS (a) 

na tc (T)a so )+ — {1 + fi P )a tc {T) a, (1) 



dT dn V dT ' au J da 

where dn{T)/dT is the thermo-optic coefficient of silica, a so = (An/n) / (AL stla i n / L$) = —0.206 is the stress-optic 
coefficient of silica [25], and [ip = —0.168 is Poisson's ratio for silica |25j . 

In the literature, the thermo-optic coefficient of silica is well characterized for room temperature. However, few 
studies have been carried out for temperatures exceeding ^800 K. To our knowledge, the widest temperature range 
(300 K — 1570 K) has been characterized using the thermally induced wavelength shift of a FBG-mirror stop-band 
AAb [26 . From the latter, the thermo-optic coefficient has been extracted for a wavelength of ^1570 nm by using 



1 dA B 
A B dT 



(1 - a so )a tc (T) + n —r^r- 



(2) 



where Ab is the central wavelength of the FBG stop-band [2Z| ■ Based on this data, we applied an established model 
that predicts the dispersion of the thermo-optic coefficient [55] to compute its value at the probe wavelength of 
Ao = 852 nm used in our experiment. By fitting the temperature-dependent refractive index change with a 2nd-order 
polynomial function (coefficient of determination: R 2 = 1 — 6 • 10 -6 ) and by differentiating the latter with respect to 
the temperature, we obtain the temperature dependent thermo-optic coefficient: 

= 9.627(9) ■ 10~ 6 + 7.74(1) ■ 10~ 9 (T - 299 K) (3) 

We compared the prediction of our fit with values extracted from several literature sources that use various methods 
|27H30| and find a relative error of the predicted refractive change index change with respect to the literature values of 
A(dn/dT) / (dn/dT) < 10 %. We note that the data from all cited works also qualitatively agrees with the temperature 
dependence predicted by Eq. [3] We also note that, to our knowledge, the thermal expansion coefficient is known up 
to 1000 K [24 . For higher temperatures we therefore employ a linear extrapolation of the two highest temperature 
values measured in |24j . 

The resulting thermal change in the effective refractive index is plotted as function of temperature in Fig. [4j In all 
the above considerations, the effect of the fiber core has been neglected. This is justified because the core constitutes 
a negligible fraction (< 0.5 %) of the fiber material which contributes to the guiding of the light in the nanofiber 
section and, in addition, has similar optical properties as the cladding material. 



The TOF radius profile 



The tapered optical fibers used in this experiment are produced from standard single-mode optical fibers with a 
cladding diameter of 125 /im. For this purpose, we employ a custom-built computer-controlled fiber pulling rig which 
allows us to fabricate TOFs with a predetermined radius profile. We chose a TOF radius profile which includes a 
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FIG. 4. Temperature dependence of the effective thermo-optic coefficient and its contributions for a fiber diameter of 500 nm. 
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FIG. 5. Simulated radius profile of the tapered section of TOF-resonator #1. 



nanofiber waist with a constant subwavelength diameter that is enclosed symmetrically by two fiber tapers as shown 
in Fig. [5} In a standard optical fiber, the light is guided by the core-cladding interface whereas in the nanofiber section 
the core has a negligible effect and the light is guided the cladding- vacuum interface. In these subwavelength optical 
nanofibers, a significant fraction of the optical power propagates outside of the fiber in the form of an evanescent wave 
that extends several hundred nanometers into the vacuum. A standard fiber exhibits a small refractive index difference 
between the core and the cladding and the fundamental fiber mode can be approximated by the weakly guided LP01- 
mode. To reach a high transmission through the TOF, this mode has to be adiabatically transformed by the taper 
transition to match the fundamental HEn-mode of the strongly guiding nanofiber [IS]. This is achieved by a taper 
profile that is composed of three linear sections followed by an exponential decrease adjacent to the waist (see Fig. [s]for 
TOF-resonator #1). In the first and third linear section of the taper, the fiber diameter is reduced with comparably 
large angles, Ox and 3 , whereas in the second linear section, the transition from a core- to cladding-guided mode 
takes place and the angle ©2 has to be chosen small enough to ensure an adiabatic transformation. 

The TOFs used in the experiment are produced from Fibercore SM800 single mode optical fibers, using the pa- 
rameters given in Tab. [I] where t\ and ri are the cladding radii at the transitions between the linear taper sections. 
Previous measurements have shown that the actual TOF radius profile agrees with the simulated radius profile within 



sample 


0i 


0i 


0i 


-^waist 


TOF-resonator #1 


5 mrad 


2 mrad 


4 mrad 


10 mm 


TOF-resonator #2 


5 mrad 


2 mrad 


4 mrad 


5 mm 



TABLE I. Parameters of the two TOF radius profiles of the two samples used in this experiment, 
the measurement accuracy of ~10 % [191 EU| . 
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Thermodynamical model 

We model the system including heat transport via radiation, heat conduction in silica, and heat diffusion through 
the surrounding gas. The assumption that no temperature gradients appear in the surrounding background gas is 
valid for mean free path lengths that are much larger than the distance from the fiber to the heat sink, such as, e.g., 
the fiber mount or the vacuum chamber wall. We estimate the pressure up to which this condition is fulfilled by 
assuming a typical distance to the heat sink in our system of 10 cm and thereby obtain p < 10~ 3 mbar. We note, 
however, that the transition to the regime in which the background gas exhibits a temperature gradient can take 
place in a pressure range that covers several orders of magnitude. For this reason, we also performed an experimental 
check of the validity of our assumption, see Fig. [2]i in the article. The system can then be described by the following 
differential equation in units of power |31] : 



c p pdtTdV = V(AVT) dV- 



k b f 2 



8nM N2 T 



P (T-T )dS-dH T 



adiativc 



dF\ 



heating 



(4) 



where dPhoating is a heat source, the remaining quantities are given in Tab. ^ and the dependencies on space and 
time are omitted for clarity. The radiated heat di? ra( jiative can be either calculated using Planck's law in combination 



symbol 


description 


value from To to the highest 
temperature covered in Ref. 


temperature range 
convered in Ref. 


Ref. 


T 


room temperature 


(294 ±0.5) K 






A 


heat conductivity 


(1.3-2.1) W/(K-m) 


(0.1 - 1650) K 


|32| 


c p 


specific heat (silica) 


(700 - 1500) J/(kg-K) 


(0 - 2000) K 


L33J 


P 


density (silica) 


2200 kg/m a 




|32| 


f 


degrees of freedom of N2 


6 






M N2 


molecular mass of N2 


4.653 • 10^ 2b kg 








Boltzmann's constant 


1.381 ■ 10"-23 J/K 







TABLE II. Model parameters and their values with references, 
with the spectral emissivity of a silica— vacuum interface (see below) or using the predictions from FED calculations 

using. 

We assume the fiber is heated via surface absorption so that the local heating power per unit length is given by the 
local surface intensity I(a(x)) of the fiber mode multiplied by the fiber circumference d x Pheat(x) — 2ir a(x) I(a(x)). 
Due to the small diameter of the nanofiber, we assume that it exhibits a constant temperature across its diameter. 
Moreover, the fiber exchanges a negligible amount of heat with the vacuum chamber which, therefore, stays at room 
temperature during the experiment. This allows us to reduce the equation to one dimension in space along the fiber 
axis (x-axis). After differentiating the resulting equation with respect to x, we then find the following differential 
equation in units of power per length : 



na 2 c r 



pd t T = 7ra 2 d x (Xd x T)- 



8irM N2 T 



p (T — To) 2tTQ — <9 x iT a diativc + ^Pheating 



(5) 



where <9-ff r adiative describes the radiated power per unit length. This equation is solved numerically for a given radius 
profile using the method of lines implemented in Mathematica with a spatial discretization of 0.1 mm along the fiber 
axis for a TOF section that extends at least 1 cm beyond either side of the waist: L s ; m = -L W aist + 2 cm. At these 
points, the fiber diameter has increased to ~10 /*m and is thus a factor of 20 larger than the waist diameter. As a 
consequence, the evanescent field and hence the heating power vanishes and the heat capacity per unit length increases 
by a factor of 400. The boundary conditions are chosen such that the fiber is thermalized to ambient temperature 
at the beginning of the simulation T(t = 0, x) = To and no heat flow occurs at the boundary of the simulation 
cell: d x T(t,±L slm /2) = 0. By solving Eq. [5J we obtain a temperature profile as a function of time and position 
along the fiber, T(t,x). We then use the latter to compute the optical path-length for discrete time-steps via the 
thermally induced effective refractive index change (see supplementary informations) and compare the results from 
this simulation to the data. 

As discussed above, the largest uncertainties that enter into the simulation are the fiber radius and the complex 
refractive index of fused silica. Therefore, we solve the Eq. [5] for the maximum and minimum values of the complex 
refractive index values for type III and IV fused silica available in literature [33] with the radius profile scaled by 90 % 
and 110 %. 
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Complex refractive index of fused silica 



The complex refractive index, h = n + ik, is the sole material property that enters the calculation of the spectral 
emissive power of thermal radiation. We take the uncertainty of the literature values into account by using all four 
combinations of the maximum and minimum of each n (refractive index) and k (extinction coefficient), respectively, 
in order to compute the total radiated power. We then use the maximum and minimum value of the latter for our 
further calculations. 

Optical fibers are made of high-quality synthetic fused silica, which is an amorphous glass, with a core with typically 
a few weight percent of dopants such as, e.g., germanium. Four types of fused silica are commercially available |34J: 
Type I and Type II are produced from crystal quartz by electric melting and from quartz powder by flame fusion, 
respectively, whereas Type III and Type IV glasses are synthetic glasses produced from SiCU in a H2-O2 flame and in 
a water-vapor free plasma flame respectively. The synthetically produced glasses feature the lowest impurity content 
and, thereby, exhibit a higher transmission than the Type I and Type II glasses, in particular in the near-infrared 
and in the ultraviolet spectral ranges. 

Since our TOF fabrication employs a H2-O2 flame for the tapering process, we consider the glass in the TOF to 
be a mixture of Type III and Type IV silica. We therefore extract the minimum and maximum values of both n and 
k, respectively, that are available for these two types of glass for our calculations. For wavelength regions in which 
we did not encounter literature values for these particular types of glass, we also use data for other types of fused 
silica. The resulting minimum and maximum contributions of the complex refractive index that we employ for our 
calculations were derived from [3J] and are shown in Fig. [6] as function of the optical wavelength. 




We note that the values given for the extinction coefficient in the visible spectrum are at the lower limit of the 
parameter range that is experimentally accessible and therefore can be considered as maximum values. Since no 
further data is available and heat radiation in the visible spectrum contributes only marginally to the total heat 
transfer rate, we do not consider smaller values. 



Emissivity of a silica vacuum interface 



We now sketch the steps that are necessary for computing the emissivity of a polished silica-vacuum interface. A 
comprehensive treatment of thermal radiation of surfaces can be found in the literature [2J [35]. The total thermally 
radiated power of an object can be calculated by using the Stefan-Boltzmann equation H(T) = e(T) cr-g. T A A, where A 
is the surface area of the object, ctb is the Stefan-Boltzmann constant, and e(T) is the total hemispherical emissivity. 
The latter can be derived using Kirchhoff's law, which states that the emissivity of an interface is equal to its 
absorptivity. For a semi-infinite silica"vacuum interface, all light that enter the silica will be absorbed. For this reason, 
the directional spectral emissivity is given by e(u, f2) = 1 — R(v : fi), where v is the optical frequency, ft = {9, 0} stands 
for the polar and the azimuthal angle, and R(is, f2) is the directional spectral reflectivity. By using Planck's law for 
the spectral emissive power, 

^ T ^¥ ex P (W)-I (6) 
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the total hemispherical emissivity can then be computed: 

J dn J duW(v,T)e(y,n) cos(0) 



e(T) 



<JbT 4 



(7) 



where dfl = sm(6)d6 dtp. The directional spectral reflectivity can be derived from Frensel's formula for absorbing 
media by averaging over both polarizations (36]: 



Pica. 



hi{v) cos(0i) — h-z{v) cos(0 2 ) 



hi(v) cos(0i) + n 2 {v) cos(0 2 ) 



hi{v) cos(0i) — rii(i/) 008(02) 



cos(0i) + n 2 (z/) cos(0 2 ) 



(8) 



where hi and 0i are the complex refractive index and the angle of incidence in vacuum, respectively, ri 2 and 2 are the 
complex refractive index and the refracted angle in silica, respectively. The angle 2 is obtained from the generalized 
form of Snell's law for absorbing media |37| : 



7 = i?e(ni) 2 - Jm(ni) 2 + sin(0i) 2 , k = (i?e(ni) 2 - Jm(ni) 2 ) sin(0i) 2 - i?e(ni) 2 Tro(ni) 2 

'sin(0i) N 



arcsm 



with 



1 ± y/l - 4k/t 



(9) 
(10) 



The total hemispherical emissivity is plotted in Fig.[7]as a function of temperature for all combinations of the maximum 
and minimum values of the real and imaginary part of the refractive index. As explained in the main text of the 
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FIG. 7. Emissivity of a polished silica-vacuum interface as function of temperature. 

manuscript, only the minimum and maximum emissivity traces (green and yellow lines, respectively) are used for 
further calculations. 
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